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In this paper we analyze the energy splitting in a generic one-dimensional 
double-well potential. By using the WKB quantization we deduce an analyt- 
ical formula for the energy splitting which is the usual Landau (1973) formula 
with additional quantum corrections. The formula is based on a linear ap- 
proximation of the potential near the turning points. We study the validity 
of this formula for the double square well potential and then we discuss the 
case of a parabolic double-well potential. 

Let us consider a one-dimensional system with Hamiltonian 

*=£ + V<*). (1) 

where V(— x) = V(x) is a symmetric double-well potential. The stationary 
Schrodinger equation of the system reads 

Hi>(x) = ( - |^ + V{x)Mx) = EiP(x). (2) 

The Sturm-Liouville theorem (see, for example, Fliigge, 1971) ensures that 
for one-dimensional systems there are no degeneracies in the spectrum. Let 
■01 and ip2 be two eigenfunctions of the Schrodinger equation 

Hipx = Exty\ and Hip 2 = E 2 vp 2 , (3) 

such that ipi(-x) = ipi(x) and ip 2 (— x) = —ip 2 (x) and E\ ~ E 2 . To calculate 
the splitting AE = E 2 — Ei, we multiply the first equation by ip 2 and the sec- 
ond by ip\ and then we subtract the two resulting equations. By integrating 
from to oo we find 

AE = ft 2 ^i(0)^ 2 (0)-^i(0)^ 2 (0) _ (4) 
2m J °° ^i(x)^2(x)dx 
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We write the eigenfunctions ip\ and ip2 in terms of the right-localized function 

Mx) = ^(Mx)+Mx)) ■ (5) 

It is easy to show that E =< ip \H\ipo >= \{E 1 + E 2 ). Then, with the 
approximation 

/ ipi(x)ip 2 (x)dx ~ - / vpl(x)dx = - , (6) 
Jo 2 Jo 2, 

we get 

AE = _^ (0)^o(0) , (7) 
m 

which is the starting formula to calculate the energy splitting. One should 
observe that this quantity is always positive, because the tail of the right 
localized eigenfunction i/jq(x) at x = has the same sign for ^(O) and its 
derivative ip' (0). Another way to see this is to realize that because of Sturm- 
Liouville theorem there are no degeneracies and that implies that all pairs 
of almost degenerate states, from the ground state up, are grouped by odd 
state above the even state. 

To determine the function ^ we perform a WKB expansion of the Schrodinger 
equation. We observe that a generic eigenfunction tp of the Schrodinger equa- 
tion can always be written as 

i 

ip(x) = exp (-a(x)) , (8) 

where the phase a(x) is a complex function that satisfies the differential 
equation 

a'\x) + (-K(:r) = 2m(E - V(x)) . (9) 



The WKB expansion for the phase is given by 

oo * 

= £( 7 )W*) • (io) 

fc=0 1 

Substituting (10) into (9) and comparing like powers of H gives the recursion 
relation (n > 0) (see Bender, Olaussen and Wang, 1977) 

n 

a'l = 2m(E - V(x)), ]T a' k a' n - k + <r" n _i = . (11) 

fc=0 



With the momentum p = ^2m(E — V(x)) the first three orders in the WKB 
expansion give 

p' , P" 3 P' 2 / 10 N 

-o=P, ^1 = "^. " = ^2-8^' (12) 

from which we obtain 



a^-ilnp, a 2 = ^ + ^J P —dx. (13) 

It follows that, up to the second order, the wave-function can be written as 

1>{x) = ± exp (1 jpdx- i h[^+ 1 -j P -dx\) . (14) 

In particular, if we call a and b the two turning points corresponding to the 
energy E, the right localized wave-function ip Q is given by 

<M0) = -^=exp(-i [ a \ p \dx + h[\\?- 2 \ + \ f I^H) - (15) 



p x nJx 4 f & Jx p° 

for < rr < a (forbidden region), and 

^o(O) = ^ exp (1 /V* - <R[^ + I [ P -dx}) + 




for a < x < b (allowed region). 

It is possible to write C\ and C 2 in terms of C by imposing the uniqueness 
of the wave-function i/j at the turning points. Following Landau (1973) we 
suppose that near the turning point x = a it is possible to approximate the 
potential locally linearly by writing 

E-V(x) = F (x-a) , (17) 

with F > 0. In this case the connections at the turning point imply that 

d = Ce'* and C 2 = Ce~^ , (18) 

and the right localized function i])q{x) can be written for a < x < b as 

= ^ cos k L fix - hl ip + sL f^i + 1) ■ (19) 

To determine C we impose the normalization condition 

1= \M%)\ 2 dx, (20) 

JO 

from which we get 

C = ^[ b ^ 1/2 - ( 21 ) 

2 Ja p 

Here we have replaced the oscillatory cos 2 (.) factor in the integrand by its 
average value 1/2, because of the rapid oscillations at small value of h. The 
final formula is given by 

A-E" = AE Landau ■ AE WKB , (22) 
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where 

2ftC 2 / 1 f a . . , x 
&E Landau = exp ( — — / p \dx) (23) 

is the usual Landau (1973) formula for the energy splitting and 

A^ fl = (l + T |^i|)exp(-/_J^|^) (24) 

is the first quantum correction. We note that higher-order WKB corrections 
quickly increase in complexity (Robnik and Salasnich, 1997) but, in principle, 
they can be calculated from the equation (11). It is important to stress that 
our formula is good if the potential is sufficiently smooth so that the linear 
approximation is valid near the turning points. 

As an example, we consider the double square well potential. In this case 
the linear approximation of the potential near the turning point is not valid. 
The potential is given by 

V(x) = V for \x\ < a 

V(x) =0 for a < x < b (25) 
V(x) = oo for \x\ > b 

For this potential we have p'(x) = p"(x) = for — a < x < a and the 
corrections to the Landau (1973) formula are zero. A naive application of 
the splitting formula gives 

2h\[E i 2a 



This formula is not correct. In fact, by using the exact wave-function 

^o(O) = D exp ( — y2m(V - E)) , (27) 
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for < x < a (forbidden region), and 

^o(O) = Aexp (^2m(E - V j) + ficxp ( - ^2m(E - V j) , (28) 

for a < x < b (allowed region), and by imposing the exact matching and 
normalization conditions (Fliigge, 1971) we find 



D 



'Vo- 


E , 


E ' 


iV - 


E , 



A = j {1 ~ i^e~ ] exp \ny 2m{Vo - E) - r E > ' (29) 

D 



B = J {1 + 1^E~ ] 6XP V^V 2m ( V ° ~ E ) + J^ E ) » ( 30 ) 



and 



Then, by applying the equation (7), we obtain: 



AhEj2m(V -E) , 2a / N 

This is the exact energy splitting for the double square well potential. It 
differs by a factor 4y" ' E(Vq — E)/V from the WKB result based on the con- 
nection formulae (18) which are not justified in the present case. 

To conclude we study a potential where our WKB splitting formula can 
be applied. We consider the parabolic double-well potential given by 

V(x) = V - x 2 for \x\ < 

V(x) =0 for y% < x < b (33) 
V(x) = oo for \x\ > b 

The turning points are at x — ±a, where a = y/Vo — E. We observe that 

fb dx fWo 1 rb i 

/ — = / dx + / dx , (34) 

J a p JVW^e y j2m(E-V + x 2 ) JWo V2mE 
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and 



p"(0), 1 



(35) 



p(0) 3 1 2m(V - E) 2 ' 



It follows that the energy splitting up to the second order in the WKB 
corrections of H reads 



In this work we have taken the first step towards a systematic improvement 
of the Landau formula (Landau and Lifshitz 1973), which is the semiclassical 
leading order energy level splitting formula for pairs of almost degenerate 
levels in double well potentials. We have calculated explicitly the quantum 
corrections up to the second order, given in equations (21-24). Our approach 
is based on the usual WKB expansion in 1-dim potentials, embodied in for- 
mulae (8-11). Thus the calculation of higher corrections can in principle be 
continued by the same method, although the structure of higher terms in- 
creases in complexity very quickly. We have also shown what happens in 
cases where the assumption implicit in the Landau formula (namely the lin- 
earity of the potential around the turning points) is not satisfied: We get 
a different result even in the leading semiclassical order, and this has been 
shown for the double square well potential in equations (26) and (32). 

It is our goal to work out a more direct WKB approach to the solution of 
the same problem, by the contour integration technique, based on requiring 




(36) 



where 




(37) 
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the single valuedness of the eigenfunction, as has been done in the case of 
a single well potential in (Robnik and Salasnich 1997). This is our future 
project. 



Acknowledgments 

This work was supported by the Ministry of Science and Technology of the 
Republic of Slovenia, and by the Rector's Fund of the University of Maribor. 



10 



References 

Bender C M, Olaussen K and Wang P S 1977 Phys. Rev. D 16 1740 

Fliigge S 1971 Practical Quantum Mechanics I (Berlin: Springer) 

Gutzwiller M C 1990 Chaos in Classical and Quantum Mechanics (New York: 
Springer) 

Landau L D and Lifshitz E M 1973 Nonrelativistic Quantum Mechanics 
(Moscow: Nauka) 

Maslov V P and Fedoriuk M V 1981 Semi- Classical Approximations in Quan- 
tum Mechanics (Boston: Reidel Publishing Company) 

Robnik M and Salasnich L 1997 J Phys A: Math. Gen. 30 1711 

Robnik M and Salasnich L 1997 J Phys A: Math. Gen. 30 1719 

Voros A 1983 Ann. Inst. H. Poincare A 39 211 

11 



